Dominos and the Gaussian free field. 
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Abstract 

We define a scaling limit of the height function on the domino tiling model (dimer model) on simply- 
connected regions in Z 2 and show that it is the "massless free field" , a Gaussian process with independent 
coefficients when expanded in the eigenbasis of the Laplacian. 

Resume 

On definit une "limite d'echelle" pour la fonction d'hauteur dans le modele des dimeres dans Z 2 . Nous 
montrons que la limite est un "champ libre gaussien", un processus stochastique gaussien dont les coef- 
ficients, dans la base des fonctions propres du laplacien, sont independentes. 

AMS classification: primary 82B20, secondary 60G15 
keywords: dimer model, domino tiling, Gaussian free field. 

1 Introduction 

A domino tiling of a polyomino P in Z 2 is a tiling of P with 2x1 and 1x2 rectangles. For a polyomino 
P let fi = jJ,(P) denote the uniform measure on the set of all domino tilings of P. 

Let U C R 2 be a Jordan domain with smooth boundary. We study uniform random domino tilings of 
polyominos P € in eZ 2 which approximate U (and using dominos which are 2e x e and e x 2e rectangles). 

A domino tiling of a polyomino P t in eZ can be thought of as a random map from eZ 2 n P £ to Z m 
the following way. Let V e = eZ 2 n P e be the set of lattice points in the polyomino P e . Let h: V E — > Z 
be a function which has the property that around every lattice square of P e the four values of h are 4 
consecutive integers ho, ho + 1, ho + 2, ho + 3, with the values on any two adjacent boundary vertices of 
P £ differing by 1. The set of such functions h (up to additive constants and a global sign change) is in 
bijection with the set of domino tilings of P t : dominos cross exactly those edges whose /i-difference is 
3. The function h associated to a tiling is called its height function fl3|| . See Figure |l[ Note that the 
height function takes values in Z, not in eZ. 

Our aim is to prove that in the limit as e — > the height function on a random tiling of P e tends 
to a random (generalized) function which has a succinct description in terms of the eigenbasis of the 
Laplacian operator on U. 

Theorem 1.1 Let U be a Jordan domain with smooth boundary in R 2 . For each e > sufficiently small 
let P e be a Temperleyan polyomino approximating U as described below. Let h e be the height of a random 
domino tiling of P t and h t be its mean value. Then as e tends to 0, h t — fe e tends weakly in distribution 
to 4/y / 7r times the "massless free field" F on U, in the sense that for any smooth function (j> on U , the 
random variable X^gy <j>{x)(h e {x) — h t (x)) tends in distribution to f v (f)F. 
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Figure 1: Height function of a domino tiling. 



For the definition of Temperleyan polyominos see below. The massless free field F on U is a 
random variable taking values in the space of distributions ^\ which are continuous linear functionals on 
the space of C 1 functions on U (with a C 1 -norm). For background on the massless free field see [juj]. It 
can be defined as follows: let {/i}j>i be an L 2 -orthonormal eigenbasis for the Laplacian A = + 
on U with Dirichlet boundary conditions (that is, /, s on dU). Let A; be the eigenvalue of Then 

where the Ci are i.i.d. Gaussian random variables of mean and variance 1. Here this expression is 
interpreted as the generalized function F satisfying, for any C 1 function (f>, 

series which converges almost surely. The expression (Q) does not define a function since the series 

diverges almost everywhere. 

Remarks. 

1. The above theorem describes the limiting value of h e — h e . The limiting average value h = lim/i e 
was computed in ||: it is a harmonic function whose boundary values are given by — times the angle 
of turning of the boundary tangent counterclockwise from a fixed basepoint. (Regarding the choice of 
basepoint, see the definition of "Temperleyan" polyomino below.) 



2. Theorem 1.1 has a well-known one-dimensional analog: Let X be the sets of random maps h from 
0, i, | . . . , 1 to Z satisfying h(0) = h(l) = and \h(^-) - h(±)\ = 1. A random element of X, when 
divided by ^Jn, converges to a random function known as the "Brownian bridge" [B. In the eigenbasis 



1 Henceforth we will refer to these objects as "generalized functions" to avoid confusion. 



2 



of the one- dimensional Laplacian -^i the coefficients of the Brownian brid ge a re again independent 



Gaussians. One difference between the one-dimensional case and Theorem L.l, however, is that the 
height function h of Theorem 0] i s unnormalized. It is therefore all the more surprising that the integer- 
valued function h of Theorem p_.l| converges to a continuous- valued object. 

3. An important open problem is to compute the distribution of the height function on a non-simply 
connected domain, even an annulus. In particular for an annulus the distribution of the height difference 
between the two boundary components (in the limit e — > 0) is unknown although it was shown in || to 
depend only on the conformal modulus of the annulus. 

4. Temperley flit] gave a bijection between the u nifo rm spanning tree process on subgraphs of Z 2 



and domino tilings. The function h of Theorem 1.1 corresponds under this bijection to the "winding 
number" of the branches of a spanning tree [^|, as first conjectured by I. Benjamini. As it is an open 
question to show that a scaling limit exists for the uniform spanning tree process [[l], 0], one might hope 
that the reconstruction of the tree from its winding numbers, which is possible for e > 0, also works in 
the limit e = 0. So far th is re mains an open problem. 
5. The result of Theorem 



1.1 



depends strongly on the choice of boundary conditions for the approximating 
polyominos P t . For even slight generalizations of these boundary conditions our methods will not work: 
see jj| for a discussion of this issue. 

6. When the region U is a rectangle U = [0, a] x [0,6], the orthonormal eigenvectors of A with Dirichlet 
boundary conditions are —r sin sin 2^p, where j, k are positive integers. So in this case the massless 
free field has independent Fourier coefficients. 

7. Most of the work to prove Theorem [^] was done in where we proved Proposition 2.2, below. 

If we consider the massless free field F to be a continuous linear functional on the space of smooth 
2-forms on U (rather than on the space of smooth functions on U) then F is conformally invariant, 
in the following sense. 

Proposition 1.2 Let to be a smooth 2-form on U and let f: V — > U be a conformal bijection. Let Fu,Fv 
be the massless free fields on U and V respectively. Let X = J (J Fu{z)u>(z) and Y — J v Fv(z)f*cj(z), 
where f*u> is the pullback of ui to V . Then the random variables X and Y are equal in distribution. 

For the proof see section 0. 



2 Background and preliminaries 

2.1 Temperleyan polyominos and approximation 

Define the (i, j)-lattice square in 1? to be the lattice square whose lower left corner is A lattice 

square is said to be even if the coordinates of its lower left corner are even. A polyomino is a union of 
lattice squares which is bounded by a simple closed lattice curve. A polyomino is even if all of its corner 
squares are even, where by corner squares we mean those lattice squares adjacent to the corners and 
containing the interior angle bisector at the corner. In particular note that an edge of an even polyomino 
P' has odd length if its two extremities are both concave or both convex corners; if the extremities consist 
of one concave and one convex corner the edge length is even. Let P be a polyomino obtained from an 
even polyomino P' by removing one lattice square 6 adjacent to its boundary, where 6 is of the same 
parity as the corners of P' . Such a polyomino is called Temperleyan, and the removed square 6 is called 
its root. In Figure [l], the polyomino is Temperleyan with root the lower left (removed) square. 

All Temperleyan polyominos have domino tilings [^, section 7]. The term Temperleyan comes from 
the bijection due to Temperley between the set of spanning trees of a rectangle in 1? and the set of 
domino tilings of a rectangular region with a corner removed Jnj. This bijection was generalized in || 
and further in |7f. 
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Let U be a smooth Jordan domain with a marked point b G dU. For each e > let P e be a 
Temperleyan polyomino in el? approximating U as follows. The boundary of P € lies within O(e) of dU, 
and the counterclockwise boundary path of P e points locally into the same half-space as the (directed) 
tangent to dU which it is near. Furthermore the root b e of P e should be within O(e) of b. 

2.2 Green's functions 

Let U be a Jordan domain with basepoint b £ dU. The Green's function with Dirichlet boundary 
conditions, or simply Dirichlet Green's function, 30(21,22), is defined to be the unique function (of 22) 
satisfying Ago{zi,z 2 ) = 5 zi (z 2 ) (the Dirac delta), and which is zero when 22 G dU, where the Laplacian 
is with respect to the second variable. This function is well-defined and when 22 is near 21 has the form 

30(21,22) = ^rlog|2 2 - Zl\ +0(1). 

The Dirichlet Green's function has the following simple expression in the basis of eigenfunctions of 
the Laplacian on U : 

Lemma 2.1 

g D {zi,z 2 ) = }^ 'r-± — . 

i>l Al 

Proof. Since the eigenbasis {fi} of A is an orthonormal basis for L 2 (U), it suffices to show that for 
each i, (fi(z 2 ), Qd(z\, 22)} = Si ^ ) . But 

(fi{z 2 ),g D (zi,z 2 )) = ^-(Xifi(z 2 ),g D (zi,z 2 )) 

Ai 

= -T-{&-fi{Z2),9D{zi,Z 2 )) 

= ^-{fi(z 2 ),Ag D (z u z 2 )) 
= ^-(fi(z 2 ),S zl {z 2 )) 

Ai 
Ai 

□ 

We will also need to define the Neumann Green's function pjv. On a bounded domain, the Green's 
function with Neumann boundary conditions does not exist. However we can define (see below) the 
function which is the "difference of two Neumann Green's functions" : for points z\ , z[ G U define 
the difference of Neumann Green's functions <?jv(zi, 22) — gN(z[, 22) to be the function of 22 satisfying 
A(pjv(zi, 22) — gN(z[, 22)) = 8 Z1 (22) — <S Z / (22), where the Laplacian is with respect to the second variable, 
and which satisfies -^{gN(zi, 22) — gN(z[, z 2 )) = 0, that is, the derivative normal to the boundary with 
respect to the second variable is zero. This function is well-defined up to an additive constant, and we 
set the constant so that g(zi, b) — g(z[, b) = 0. 

Near 21 we have g N (zi,z 2 ) - g N (z' 1 ,z 2 ) = log \z 2 - 21 1 + ci(zi, z[) + 0(z 2 - 21). Let g N (zi,z 2 ) - 
cIn(z'i, 2 2 ) be the harmonic conjugate (with respect to the second variable) of gjv(^i, z 2 ) — gN(z'i, z 2 ). This 
function is multiply- valued, increasing by 1 when 22 turns counterclockwise around 21 and by —1 when 
22 turns counterclockwise around z[. The function §n(zi,z 2 ) — 3^(2^,22) := gN(zi,z 2 ) — gN(z[,z 2 ) + 
i(g N (zi, 2 2 ) — gN{z[,z 2 )) is analytic in z 2 (except at 21 and z[) and is the analytic Neumann Green's 
function. It is also multiply-valued. 

We can define the exterior derivative of §n with respect to the first variable as d§N(zi,z 2 ) = 
Fq(zi, z 2 )dx\ + iF 1 (x,z 2 )dyi, where z\ = X\ + iyi and F ,Fi are defined by taking limits (for <5 real) 
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Fo( Zl ,z 2 ) = lim^o ™(*i+s.*2)-§»(zi,'2) and iFl ( Zl ,z 2 ) = lim^o M<£l+i^2t±iiSl^ , 
tions Fq,Fi are well-defined (single- valued) and vanish at 22 = 6. 

As examples of these functions, on the upper half-plane H with 6 = 00 we have 



g D (zi,z 2 ) = 7~; lo g 



22 - 21 



These func- 



(2) 



and 



g N (zi,z 2 ) - g N {z[,z 2 ) 



1 , («2 - 2l)(Z2 - Zl) 

— log — ==- . 

27T & ( Z2 - Z[)(Z 2 - Z[) 



(3) 



Note that when Z2 £ E the imaginary part of (^) is constant (in fact vanishes); this implies that the 
real part has Neumann boundary conditions. For a more general Jordan domain V , let / be a Riemann 
map from V to the upper half-plane sending b (the base point of V) to 00. Then the Dirichlet Green's 
function on V is ffK (21,22) = <to(/(0i), f(zi)), and the analytic Neumann Green's function is defined 
similarly g%(zi,Z2) — 3^(2^,22) = <?Jv(/(2i), f{z2)) — (f(z'i), f( z i))- One can in fact take this to be 
the definition of the Green's functions on V. 



2.3 Moment formula 

For a region U with basepoint b £ dll, define the functions F+(zx, 22) and F-(zi, 22) by 

-4dg N (zi, 22) = F+(zi, z 2 )dzi + F-{zi, z 2 )dzT, 



(4) 



where d is exterior differentiation with respect to the first variable. Then F+(z\,b) = = F-(zi,b) (in 
terms of the functions Fo,F\ of the previous section we have F± = — 2(Fb ± -Fi)). 
Let ho(x) = /i(i) — h(x). 

Proposition 2.2 (|H[) Under the hypotheses of Theorem \l . ^ , let z 1: ... ,2& 6e distinct points ofU, and 
71, . . . , 7fe disjoint paths running from the boundary of U to zi, ■ ■ . , 0ft respectively. Let h(zi) denote the 
height of a point of P t lying within 0(e) of Z\. Then 



lim E(/lo(2i) ■ • ■ ho(z k )) 



ei,...,e*e{±l} 



/ ••■ / det (F £l , £j ( 2l , 2j ))d 2 < £l) ...d4^ ) , 



(5) 



where dz^ = dzj and dz\ 1 ' = rfzj, and 



( 



if i = 3 



F + {zi,z 3 ) if ( Ei)Ej ) = (1,1) 

i ? e i ,e,(2 i ,0 J ) = ^ J 1 -^,^-) if (Ei,Ej) =(~1,1) 

F-{zj,Zj) if (ei,£i) = (1,-1) 

F+(zi, Zj ) if ( ei) £,-) = (-!, -1). 



3 Proof of Theorem 



When U is the upper half plane with basepoint at 00, the derivative of the analytic Neumann Green's 
function is (see (|3|)) 

dzi dzT 



dgn(zi, 22) = 



27r(zi — 22) ^ 2ty(zi — 22) 



Thus from 



we have F+(zi, 22) = — -, — - — 7 and F- (21, 22) = —7 — 2 > 

' ' 7r(z 2 — Zl) V x > 7r(z 2 -zi) 
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Let p,q£U. From Proposition 2.2 we have lim e _»o K(ho(p)ho(q)) 





F+ (zi,z 2 ) 
F+(z 2 ,zi) 



F- (zi, z 2 ) 

F- (z 2 ,zi) 



dzidz2- 



' 71 >72 

Plugging in for F± gives 



F-( Zl ,z 2 ) 
F-(z2,zi) 



dz\d~z~2 + 



lim E(/io(p)fto(g)) = — 4r / / 7 — rrrdzidz2 H 5- / 

E ^o v/v wy n 2 7 71 7 72 (z 2 - zi) 2 n 2 J 11 



^+(21,^2) 
F+(z2,zi) 



dz\dz 2 . 



(22 — 21) 



+ - 



1 , ,_ 4 

— —dzidz^ 

'- TV 2 



71 ^72 ( 22 - Zl ) 

g 

= —Re log 



-— — ——dzTdz2 

(22 - Zi) 2 



p-q 
p- q 



(6) 



Note that this is — ^§go(p,q) where qd is the Dirichlet Green's function on U (see (Q)). 
Now let pi, . . . ,pk be distinct points in the upper half plane U. By symmetry of ho, if k is odd the 
moment E(/io(pi) • ■ ■ ho(Pk)) is zero. We therefore assume k is even. In (H), the matrix has y entry 

o 

''<:« ,Cj (2i, 2j) = 



*Vj -2- ") 



Such a matrix has a simple determinant: 

Lemma 3.1 For k even let M be the k x k matrix M = (rriij) with mu — and rriij = i iw/ien 

det(M) = ^ — , (7) 

where the sum is over all (k — 1)!! possible pairings {{a(l), cr(2)}, . . . , {cr(fc — 1), cr(fc)}} 0/ {1, . . . , k}. 

This lemma also appears in ^j. 

Proof. The proof is by induction on k. The formula clearly holds when k = 2. For k > 2, the 
determinant is a rational function of Xi with a double pole at xi = x 2 ; we can write 



det(M) 



C-2 



+ 



C-l 



(xi - x 2 ) 2 {xi - x 2 ) 



CO + 0(xi — x 2 ) 



The coefficient c_i is zero since the determinant is even under the exchange of x\ and x 2 (exchange the 
first two rows and exchange the first two columns). The coefficient c_2 is the determinant of M12, the 
matrix obtained from M by deleting the first two rows and columns. Therefore the right and left-hand 
sides of (Q) both represent rational functions (in each variable) with the same poles and residues; hence 
they differ by a constant. This constant is zero by homogeneity. □ 



Combining the lemma with Proposition gives the following. 

Proposition 3.2 Let U be a Jordan domain with smooth boundary. Let pi, . . . ,Pk £ U (with k even) 
be distinct points. We have 



JimE(fto(pi) ■ • ■ ho(p k )) = (~~~) 3D(p CT (i),P CT (2))---ffD(p CT (fc-i),f> CT ( fc )) 

pairings a 



k/2 
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3.1 



Proof. When U is the upper half plane this follows by combining Proposition 2.2 with Lemma 
and the calculation (^). For arbitrary U, equation Q shows that E(/io(f>i)/io(P2)) = — ^§gr>(pi,P2) 
(where g% is the Dirichlet Green's function on U) by conformal invariance of the height moments and of 
go- This completes the proof. □ 



The proof of Theorem 1.1 is completed as follows. Let f ni , . . . , f„ k be (not necessarily distinct) 



eigenvectors of A with Dirichlet boundary conditions. Let C„) be the real-valued random variable 

Cnj = e 2 52 x ev ho(x)f nj (x) , where the sum is over the vertices V e of P t , and f n . (x) is f nj evaluated at 
the vertex x. We have 
limE((?W---C^) = 



lim •■■ e 2 f ni (x 1 )---e 2 f nk (x k )E(h (x 1 )---h (x k )) 

e— >0 z — J z — J 



fc/2 



j "'/ /™l( X 0'--/^( :r fc)5Z 5l5 ( :E<T ( 1 )' a::<T ( 2 ))''' 5D ^' T (' i: - 1 )' :E<T (' i: )) 

16\ fc/2 ^ / [ f I \ t I \ /m 1 (a CT (l))/ mi (a: CT (2)) /m fc/2 {Xg(k-l))fm k/2 (x a ( k )) 

*/ a J U JU mi ,..., mk/2 Am l Am M2 

16^ fc,/2 <5<t(1),<t(2) 5o-(fc-l),<7(fc) 



E 



(-A CT (1)) (-ACTffc.!)) 



By Wick's theorem [|10| , these are exactly the moments for a set of independent Gaussians of mean zero 
and variances — Now to conclude we invoke the following standard probability lemma : 

Lemma 3.3 ([Q]) A sequence of (multidimensional) random variables whose moments converge to the 
moments of a Gaussian, converges itself to a Gaussian. 

This completes the proof. 



4 Proof of Proposition 



Since X and Y are Gaussians (each being the sum of Gaussians), and have mean 0, it suffices to compute 
their variances. But 

Epf 2 ) = f ( lu{z 1 )lo(z2)E(F(z 1 )F{z 2 )) 
Ju Ju 

L0{z 1 )u](z 2 )g r \{z\, Z 2 ) 

/*w(yi)/*w(y2)5fl(/(yi), f{y2)) 

V JV 

/*w(j/i)/*w({/ 2 )flc(j/i,3/a) 



IV Jv 
= E(Y 2 ) 

where we used the conformal invariance of the Green's function, go(f(yi), /(j/a)) = <?z>(j/i , 2/2). This 
completes the proof. 
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